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Abstract 

We present a derivation of a formula that gives dynamics of an integrable cellular automaton 
associated with crystal bases. This automaton is related to type D afhne Lie algebra and contains 
usual box-ball systems as a special case. The dynamics is described by means of such objects as 
carriers, particles, and antiparticles. We derive it from an analysis of a recently obtained formula 
of the combinatorial R (an intertwiner between tensor products of crystals) that was found in a 
study of geometric crystals. 

I INTRODUCTION 

The crystal basis theory [10, 11] has played an important role in studies of solvable lattice mod- 
els and integrable systems since more than a decade. In this context Ref. [9] by Kang, Kashi- 
wara and Misra has provided useful families of crystal bases associated with affine Lie algebras 
A£\B£\c£ ) ,D£\A!£ t l-i,A < £{ and D^. We call the first algebra type A and the fourth one type D. 
In Ref. [3] we have obtained a formula of the combinatorial R (an intertwiner of the crystals) associated 
with the type A algebra. From the viewpoint of integrable systems an intriguing fact in Ref. [3] is 
that this formula has been derived from a discrete soliton equation (the nonautonomous discrete KP 
equation) by a procedure known as the ultradiscretization. Since then there has been progress which 
goes beyond the type A case. In a study of geometric crystals associated with the type D algebra 
we have obtained an explicit formula of a tropical R, an intertwiner of geometric crystals [13]. The 
tropical R is a birational map between totally positive rational functions, while the combinatorial R is 
a bijective map between finite sets. Further analysis of this tropical R and the combinatorial R derived 
from it should be an important task in studies of integrable systems, since they are connected with 
discrete and ultradiscrete soliton equations of type D Lie algebra symmetry [14] . 

The purpose of this paper is to investigate a piecewise linear formula of the above mentioned 
combinatorial R for the type D crystals in Ref. [13]. Our main result is the derivation of the limit 
of the formula that leads to the particle antiparticle description of an integrable cellular automaton 
(Theorem 13). This description was recently obtained [15] by using a factorization of the combinatorial 
R into Weyl group operators, a property that had been found and proved in Ref. [7]. We emphasize 
that the result proves a non-trivial fact that the factorization of the combinatorial R can be conducted 
in two different ways, via the Weyl operator description and via the piecewise linear formula. This 
point is new even in the type A case (Theorem 3). 
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We briefly explain the background to our problem. There were studies on one dimensional cellular 
automata known as the box-ball systems [17, 18, 19, 20, 21]. It was found that dynamics in these 
automata was controlled by the combinatorial R of the type A crystals [2, 3]. Based on the families 
of crystals in Rcf. [9] integrable cellular automata associated with crystals of the other types were also 
constructed [5, 6]. A question about such generalized automata arose as to whether we could give a 
description of their dynamics as box-ball like systems. To answer this question the particle antiparticlc 
description was found [8, 15]. 

In Rcf. [15] it was found that the automata associated with crystal bases of any types of affine Lie 
algebra in Ref. [9] can be embedded into the type D case. Thus one can obtain the particle antiparticlc 
description of these automata from that of the type D case. This is the reason why we devote ourselves 
into this particular case. 

The plan of this paper is as follows. In Sec. II the automaton associated with the type A crystals is 
reviewed. The piecewise linear formula of the combinatorial R is presented. The particle description of 
the automaton is proved in terms of the piecewise linear formula. In Sec. Ill we discuss the piecewise 
linear formula of the combinatorial R of the type D crystals. Its reduction to the type A case is also 
shown. The automaton associated with the type D crystals is explained in Sec. IV. The formula of a 
factorized dynamics of an inhomogeneous automaton (Theorem 13) is reviewed. This formula is proved 
in Sec. V using the piecewise linear formula of the combinatorial R. Proofs of several lemmas are given 
in the Appendix. 



II A^j CASE 

A Combinatorial R 

We begin with type A case. Instead of A$ we adopt A^li crystals because it enables us to compare 
the results with those in the crystals. For the notation we use overlines to distinguish the symbols 
from those in the type D case, writing B for a crystal, R for the combinatorial R and so on. 
As a set the A^\ 1 crystal Bi (I is any positive integer) is given by 



Bi = I {xi,--- ,x n ) ezi 



X> = 4- (!) 



2 = 1 



The other properties of this crystal are available in Ref. [9] . In this paper we use no other property of 
Bi. We will write simply B or B for Bi with arbitrary I. 

Definition 1. Given a pair of variables x = (xi, . . . , x n ) € B, y = (j/i, . . . , y n ) S B , let R : (x, y) i— > 
(x 1 , y') be the piecewise linear map defined by x' = (x[, . . . , x' n ), y' — (y[, . . . , y' n ) where 

x'i = Vi + Pi+i - Pi, 

y'i = Xi + Pi- p i+1 . 

Here Pi is given by 

/j-i \ 



Pi = max V7t/ fe+i _i - x fe+i _i) + y j+i -i . (2) 

Kj<n \ ' / 
\fe=l / 



The indices herein involved are interpreted in modulo n. 
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Except for the notation this formula is the same one as defined in Proposition 4.1 of Rcf. [3]. The 
normalization of (2) is so chosen as the Pi to take the same expression as the formula in Theorem 5.1 
of Rcf. [9] . The property of R to intertwine the actions of Kashiwara operators in the crystal basis 
theory was essentially proved in Sec. 1 of Ref. [13]. It ensures that the (x',y') falls into B x B. 

We note that there exist other ways to present the combinatorial R (e.g. Refs. [4] and [16]) but 
which are not used in this paper. 

B Automaton 

Now we consider the automaton. In Ref. [3] the space of automaton extends infinitely towards both 
ends 

• • • x Bi^, x B h x ~B k+1 x • • • . 
For our purposes it is sufficient to consider a finite size system like 

B h x ••• xB, r (3) 

Let L{~^> J^iLi '*) be an integer. We call a particular letter for the ground state of the automaton a 
vacuum. We can use any letter in {1, . . . ,n} as a vacuum of the automaton [15]. Throughout this 
paper we adopt n as the letter for the vacuum. As in Ref. [7] we define 

Bh[n] = \ (xi, . . . ,x n ) e B L 



x n 3> x a for any a ^ n 

We write B x B' ~ B' x B for correspondence by R. Take any x <E £?l [«] • Applying R successively we 
have 

B L \ n] x (Bi ± x • • • x -Bjjy) — {Bi x x • • • x -Bz N ) x _B^[n], 

xxY ^ X' xy', (4) 

that gives the following. 

Definition 2. The time evolution operator T of the automaton is given by 

T :Y i-> X'. 

It means that we regard Y and X' in (4) as two automaton states before and after the time evolution. 
We note that the operator T actually depends on x. 



C Particle description 

There is an interpretation of the automaton that we call a particle description. It is the description 
of the box-ball systems in Refs. [17, 18, 19, 20]. Suppose we have balls with index a (1 < a < n — 1) 
that we call a-balls. For x — {x\, . . . ,x n ) G Bi i we associate a box of capacity U that has x a a-balls 
(1 < a < n — 1) in it. Then an element of Bi ± x • • • x Bi N is regarded as a one dimensional array of 
boxes of capacities li, . . . , In with these balls. For any a (1 < a < n — 1) we consider a carrier of a-balls 
that we call an a-carrier. We assume that the a-carrier has a sufficiently large capacity, so that it can 
carry arbitrary number of a-balls at a time. 
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First we suppose h = 1 for all i. We call the associated automaton basic [15]. In this case the x 
represents a box with an a-ball i£ ar = 1 for a ^ n. It represents an empty box if x n = 1. For any a 
we write j^aj for x with x a = 1. The carrier goes along the array of boxes. Then there arc four actions 
in the loading-unloading process by the a-carrier: 

1. If the carrier has at least one ball and meets an empty box, we unload a ball from the carrier 
and put it into the box. 

2. If the carrier meets a box with an a-ball, we pick up the ball and load it into the carrier. 

3. If the carrier meets a box with a &-ball (b ^ a), we do nothing. 

4. If the carrier has no ball and meets an empty box, we do nothing. 

These actions are depicted by the left four pictures in Fig. 1. For any a (1 < a < n — 1) let fC a be a 
particle motion operator that acts on the space of automaton (array of boxes) and does the actions in 
the loading-unloading process explained above. We assume that the JC a depends on x in (4) in such a 
way that the a-carrier has x a balls in it at the beginning where x a is the a-th element of the x. Then 
for the basic automaton we have [2, 3] 

T = K,\K,2 ■ ■ ■ fcn-l- 

Now we consider a not necessarily basic case which we call inhomogeneous [15]. We denote by 
V the operator that reduces the automaton into a basic one, and by Q the operator that makes a 
rearrangement of balls [15]. To explain them we first let N = 1 in (4) 

S L [n]xB,~I,xB t [n]. 

Then let 

V : Bi -> Bi x • • • x Bi 

S v ' 

I 

be the operator which sends y = (t/i, . . . , y n ) £ Bi into 

p7|x---xp7|x---x[Tjx---x[Tj. 

S v ' S v ' 

yn yi 

Its inverse V~ x can be defined only on such arrays in which the letters are arranged in decreasing order. 
Let 

Q : Bx x • • • x B 1 -> ~B X x • • • x ~B X 

" v ' S v ' 

/ / 

be the operator which packs |~nj 's into the left end. Next we consider the case TV > 1. We insert walls 

between the Bi i 's in (3) to mark their positions. Then by V we denote the operator that applies the 
above V on each B^, and by Q or V^ 1 those that applies the above Q or V~ x on each B\ x • • • x B\ 
between the walls. Now we have [1, 15] the following. 

Theorem 3. The time evolution operator of the inhomogeneous automaton is given by 

T = V- 1 QK 1 K 2 ...K, n - 1 V. (5) 

This theorem means that the time evolution of the inhomogeneous automaton can be reduced into 
that of the basic one only by inserting a simple rearrangement. 
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D Proof of the particle description 

The reduction of an inhomogeneous automaton to a basic one (Theorem 3) was first presented by 
Fukuda [1]. We gave another proof of this theorem (and its generalization to type D case) in Ref. [15]. 
Here we show still another proof of this theorem that uses the piecewise linear formula in Definition 
1. This is the proof that was inferred in Sec. II D of Ref. [3] but was not explicitly given there. Let 
Pi = linxj^^oc Pi. Then by (2) we have 

Pi = ma x . V(j/fc+*-i - x k+t -i) + yj+i-i . (6) 

From Definition 1 and (6) we obtain p n = y n and 

p % = max{ ? y t , y t - x t + p i+ i}, (7) 
x'i = minfe+i,^}. (8) 

Note that the relation (7) is a descending recursion formula for p^s on i. Let x G Bl[ti] and y G Bi 
be a pair of variables. Let K, a be the particle motion operator introduced in Sec. II. C. Then it is easy 
to see that the piS (respectively x'^s) for 1 < i < n — 1 obtained by (7) (respectively by (8) ) denote 
the number of empty boxes (respectively the number of boxes with balls with index i) in the automaton 
state K,i . . . lC n -iPy. We also see that the a-carrier finally has y' a = x a + y a — x' a balls in it. This proves 
the factorization of T in (5) for the case N = 1 in (4). The assertion of the theorem for the case N > 1 
follows immediately by repeated use of this case, where we adopt the final states of the carriers for Bi i 
in (3) as their initial states for Bi i+1 . 

Ill dw combinatorial r 

A Piecewise linear formula 

As a set the crystal Bi (I is any positive integer) is given by 



Bi = I (x ll ■ ■ ■ ,x n ,x n , ■ ■ ■ ,x{) G Z' 



x n x n =0,^2(xi + Xi) = l\ . (9) 



i=l 



The other properties of this crystal are available in a preprint version of Ref. [9] (Kyoto Univ. , RIMS- 
887, 1992) or in, e.g., Ref. [12]. In this paper we use no other property of £>;. We will write simply B 
or B' for Bi with arbitrary I. 

Definition 4. Let x — (x\, . . . ,x\) G B,y = (j/i, . . . ,y x ) G B' be a pair of variables. The involutive 
automorphisms *, a\, a n on x, y are defined by 

* :Xi <-+yi,Xi <-+ yi (1 < i < n), 
CTi :xi <-+x 1 ,yi 
o~n '■ x n <-> x n , y n <— > y n . 

For any function F = F(x, y) we denote by F a the function obtained from F by applying a — (*, a\, cr„) 
to it. For i e Bi we write £(x) for /. 
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Definition 5. Given a pair of variables x = (x\, . . . , xi) G B, y = (j/i, . . . , y 1 ) e £>', let i? : (x, y) 
(x 1 , y') be the piecewise linear map defined by x' = (x[, . . . ,x[), y' = (y[, . . . ,y[) where 



x[ 


= 2/1 


+ 


Vo 1 










x 'i 


= Vi 


+ 


Vi-! 


- ^i + Wi 


-Wi-i 


(2 < i < n- 


1), 


x n 


= y n 


+ 


Vn-1 


- C", 








x\ 


= Vi 


+ 








(1 < i < n), 




y'i 


X{ 


+ 


v t — 1 


- v;* 




(1 < i < n), 




y'i 


= X\ 


+ 


yai 
v 


- V{, 








% 


~ X% 


+ 


y %— i 


- V? + Wi 


- Wi-i 


(2<i<n- 


1), 


Vn 


-- x n 


+ 


V* , 

^n— 1 


" n ■ 









(10) 



Here Vi and Wj are given by 

V"i = max {a it j,a'i ■}, 

l<j<n— 1 

VK, = max (V5 + - y u Vi-i + V* - x,) + mm{ Xi , %), (1 < i < n - 2), 
W n _i = F n + V^». 

The functions atj — aij(x,y) and a-j = a' i j(x,y) in (11) are given by 



&i,j{x,y) = max(<5 J> _i/3 i ,y J - Xj) + 



*(v) + Efc=i+i (** - wit) for .? > *. 



= max^-i/J-,^ - y^) + l{x) + ^(y fe - x k ) + ^(j/k - ar fe ) 



(11) 

(12) 
(13) 



(14) 
(15) 



fe=i 



fe=i 



where 



i„ — y n for i n — 1, n, 

fori = n-l, /3- 

x n -y n forz = n, 



y„ — x„ for i ^ n — 1, n, 

max(y„ - 2x„, y„ - 2a;„) fori = n-l, 
y n -x n for i = n. 



The map i? is the combinatorial i? for the D„ crystals. The property of R to intertwine the actions 
of Kashiwara operators in the crystal basis theory was proved in Theorem 4.28 of Ref. [13]. It ensures 
that the (x\ y') falls into B' x B. 

Remark 6. We have changed the notation from Ref. [13] since our present formalism uses both x n 
and x n . The changes are listed in Table I. 

According to the correspondence in Table I one of the formulas in Eq. (4.66) of Ref. [13] is now 
translated into 



Vn + max(K - y n , V° n - y n ) - V° n , 
x' n = y n + max(K - y n , V° n - y n ) - V n . 



x 



In order to make them coincide with the relations in (10) we should define the V n -\ as max(V r „ — 
IJm V° n — Vn)- Actually the above definition of V n -\ is equivalent to this. In other words we have the 
following. 
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Table I: The correspondence of the notation between in this paper and in Ref. [13] for the piecewise 
linear formula of the D n ^ combinatorial R. The z in the first column denotes x,y,x' or y' . It is 
assumed that i ^ n — l,n in the last two columns. 



This paper 


Ref. [13] 


This paper 


Ref. [13] 


This paper 


Ref. [13] 




max(z„, 0) 




Vn 


K-l 


&i,j(^n— 1) 




%n 


max(— z n , 0) 






v n—1 




max (^,j>1ij) 


Z n -1 


z n -i + mm(z n , 


0) 








max(»7 iin _i,77 iin ) 


z n-l 


z n -i + min(z„, 


0) 




* o a n 


<»-! 





Lemma 7. The following relation holds: 

V n -i = max(K - y n , V° n - y n ). (16) 

Proof. We have 

a n-i.j = max(a„j - y„, - j/ n ), 

= - (a™j) CT " - Vn), 

for 1 < j < n — 1. In order to check these relations we can use max(— x n , —x n ) — max(— y n , —y n ) — 0. 
The claim of the lemma follows immediately from these relations. □ 

Remark 8. The transformation properties of the piecewise linear functions Vj and Wi under the 
automorphisms 01, *, and u n will be used afterwards, so we list them in Table II. It was quoted from 
Ref. [13] and adjusted by the correspondence in Table I. 



Table II: The transformation of the piecewise linear functions V% and Wi in Definition 5 by the auto- 
morphisms <7i, *, and er„ in Definition 4. 





V 


V l (Ki<n-l) 


V n 


Wi{l<i<n-l) 


o-i 


v 


Vi 


V n 


Wi 


* 


V 


V* 


V n 


Wi 


0~n 


V 


Vi 


" n 





Some more relations on the piecewise linear functions will be used later. We give them at the beginning 
of the Appendix. 

B Reduction to the case 

We realize that the piecewise linear formula in type D case has a rather bulky expression in contrast 
with its type A counterpart: See Definitions 1 and 5. In order to understand its structure it is worth 
trying to study some special limits of the formula. Here we consider a reduction to the type A 
case. We observe that the piecewise linear map R in Definition 5 for the crystals becomes the 
intertwiner of the A^li crystals under the reduction. 
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Theorem 9. Set 

x i = y i = for l<i< n. (17) 
Then the map R : (x,y) i— > (x',y') in Definition 5 reduces to 

x • = y t + P i+1 - Pi, 
y'i = Xi + P { - P l+1 , 
x'i = 0, 
Vi = 0, 

where Pi was defined by (2). 

This theorem follows from Lemma 10 below. 
Lemma 10. Under the specialization (17) the following relations hold 

Vi=t{x)+P u V*=£(x)+P i+1 (0<i<n), 

V° 1 =l(x) + P 2 , V^=t{x) + P n , 

W i = 2l{x) + P l +P i+1 (l<i<n-l). 

We shall give a proof of this lemma in Appendix. 

We note that the reduction from type D to type A (Theorem 9) itself can also be obtained from 
the description of the combinatorial R in Ref. [4] since the insertion algorithms for types A and D in 
Rcf. [4] coincide under the condition (17). 

IV dw automaton 

A Definition 

We now present a brief definition of the automaton using the crystals and the combinatorial R. 
For a more complete definition, see Refs. [5] and [6]. We consider a finite size system like 



B h x ••• x Bi N . 



Let L{~^> J2iLi h) DC an integer. We define 



BlM = < (xi, . . . ,x n ,x n , . . . ,xi) e B L 



x n ^ x a for any a ^ n > . (18) 



Take any x E Bl[u\. Applying the combinatorial R successively we have 

Bj\ n] x (Bi 1 x • • • x Bi N ) ~ x • • • x Bi N ) x Bl[u], 

x x Y m. X' x y', (19) 

that gives the following. 

Definition 11. The time evolution operator T of the automaton is given by 

T :Y \ — > X' . 

Remarks similar to those after Definition 2 also apply here. 
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B Particle antiparticle description 



We consider a particle antiparticle description of this automaton [8, 15]. This is a generalization of 
the particle description in Sec. Ill C. Suppose we have balls with index a and a, (1 < a < n — 1) 
that we call an a-ball and an a-ball respectively. The a-ball and the a-ball are regarded as a particle 
and an antiparticle one another. We introduce a pair annihilation process in which a pair of particle 
and antiparticle makes a bound state, and a pair creation process where the bound state breaks up 
into a pair of particle and antiparticle of another kind. To each x — (xi, . . . ,x n ,x n , . . . ,x\) <E Bi i we 
associate a box of capacity that has x a a-balls, x a a-balls (1 < a < n — 1), and x n bound states in it. 
Then any element of Bi x x • • • x Bi N can be regarded as a one dimensional array of boxes of capacities 
l\, . . . , In with the balls and the bound states. For any a (1 < a < n — 1) we introduce the notion of 
an a-carrier as in Sec. Ill C, and that of a carrier for a-balls that we call an a-carrier. Assume that 
their capacities are sufficiently large. 

First we consider a basic case, i.e. we suppose k = 1 for all i. In this case the x represents a box 
with an a-ball if x a = 1 and a box with an a-ball if x a = 1 for a ^ n. It represents an empty box if 
x n = 1 and a box with a bound state if x n — 1. We write j^aj for x with x a = 1 and write for x 
with x a = 1 . 

Remark 12. In what follows we write a also for a number with an overline as well as that without an 
overline. We interpret a — a and x-a ~x a . We call x a the a-th element of x. 

Besides the four actions in the loading-unloading process in Sec. Ill C we have three additional 
actions by the a-carrier: 

5. If the carrier has at least one ball and meets a box with an a-ball, we unload a ball from the 
carrier and make a bound state in the box. 

6. If the carrier meets a box with a bound state, we extract an a-ball from the bound state, load it 
into the carrier, and leave an a-ball in the box. 

7. If the carrier has no ball and meets a box with an a-ball, we do nothing. 

These actions are depicted by the right three pictures in Figure 1. 

For any a £ {1, . . . , n — 1} U {n — 1, . . . , 1} (see Remark 12) let JC a be a particle motion operator 
[8, 15] that acts on the space of automaton and does the actions in the loading-unloading process 
explained above. We assume that the K, a depends on x in (19) in such a way that the a-carrier has x a 
balls in it at the beginning, where x a is the a-th element of x. Then for the basic automaton we have 



Now we consider an inhomogeneous case. We define the operators V and Q as in Sec. Ill C but 
modify them to be suitable for the type D case [15]. To define them we first set N — 1 in (19) 



[8] 



•n— 1- 



B L [n] xBi-BtX B L [n]. 



(20) 



Then let 



V : B t 



B 1 x • • • x B x 
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be the operator which sends y = (j/i, . . . , y n , y n , . . . , j/J 6 Bi into 

x • • • x j"Tj x • • • x x • • • x x x • • • x |~n~| x • • • x |TJ x • • • x |TJ. (21) 

V -v ' v -V ' v v ' V V ' 

Vl Vn Vn VI 

Its inverse V~ x can be defined only on such arrays in which the letters are arranged as in (21). Let 

Q : Bx x • • • x Bx -> B 1 x • • • x B t 

" v ' S v ' 

I I 

be the operator which packs |~rej 's into the left end and p^j 's into the right end. For the case N > 1 
we generalize the definitions of these operators in the same way as in Sec. Ill C. Now we have the 
following. 

Theorem 13. The time evolution operator of the inhomogeneous automaton is given by 

T = V^IC— . . . %/C T Q/CiX; 2 . . . Kn-xV. (22) 

In Ref. [15] a proof of this theorem was given by means of the factorization of the combinatorial 
R in Ref. [7] . In the remaining part of this paper we give another proof of this theorem that uses the 
piecewise linear formula of the combinatorial R in Definition 5. 



V PROOF OF THE PARTICLE ANTIPARTICLE DESCRIP- 
TION 

A Limit of the piecewise linear formula 

We study a limit of the piecewise linear formula of the map R in Definition 5. Let F = F(x, y) be any 
function of (x,y) E B x B' . The limit we consider here is to adopt the -Bz,[n] in (18) as the B. We 
introduce the following normalized limits: 

lim*F= lim (F(x,y)-£(x)), (23) 

x n ^oo,x n —*0 

linwF= lim {F(x,y)-2£(x)). (24) 

First we consider the limit (23) of Vi . For the sake of notational simplicity we denote lim* Vi by Vi . 
For a = <Ti , a n or * we denote lim* V® by vf. Note that if a = a n or * the vf is not necessarily equal 
to the function that is obtained from Vi by applying a to it, whereas if a = a\ it is. Next we consider 
the limit (24) of Wt . We shall denote lim** Wt by to, . 

The relations in Lemma 23 in the Appendix become recursion relations in the limit (23). We let 
(x) + denote max(x,0). 

Lemma 14. For F = Vi, V* , Vq 1 or V° n the limit lim* F exists. Moreover the following relations hold 



K n = v n-l=yn-y n , (25) 

v*_ 1 =y i -Xi + max.{v*,(x i -y i ) + }, (1 < i < n - 1) (26) 

v i = max{y i -x i +Vi-i,(y i -Xi) + }, (1 < i < n - 1) (27) 

v n = y n -Vn+ max{y„ + y n _ x -x n -i + v n - 2 , (y n + y n -i - x n -i)+}- (28) 
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We consider the limit (24) of the defining relations of Wi ((12) and (13)). 

Lemma 15. The following relation holds 

w i ^v l + v* - min{v* - v*_ x + y l7 v l - Vi-i + x t } + min{x;, yj (1 < i < n - 1). (29) 

These lemmas will be used in proofs of Lemmas 20 and 21. We shall give their proofs in the 
Appendix. 

B Analysis of the particle antiparticle description 

We now consider a recursion formula satisfied by the numbers of items in the particle antiparticle 
description. For this purpose we introduce the following. 

Definition 16. For any non-negative integers A, B, C, D, E, define the piecewise linear map 

7 : (A, B, C, D, E) h-> (F, G, H, I, J) 

by 

F = mm(A, E), 

G = B + {A-E)+, 

H = min (C,B + (E-A)+), 

I = D + (C-B-{E-A)+) + , 

J = D + {B-C+(E-A)+) + . 

The identities F + G = A + B, H + I = C + D, and F + H + J = B + D + E can be checked 
easily and will be used afterwards. We give an interpretation of the map 7 in the particle antiparticle 
description that is illustrated in Figure 2. Recall the seven actions in the loading- unloading process 
by the a-carrier (Figure 1) that were explained in Sees. II C and IV B. We write act-i for the action 
with number i in the lists. Note that the a can represent an overlined number in Figure 1. In Figure 2 
the boxes with 6-balls (b 7^ a, a) have been omitted because of act-3. In the upper picture of Figure 2 
we are applying act-5 (or act-7 if E = 0), act-6, act-1 (or act-4 if B + (E — A) + = 0), and act-2 from 
left to right. In the lower picture we are applying act-1 (or act-4 if E = 0), act-2, act-5 (or act-7 if 
B + (E — A) + = 0), and act-6. 

In what follows we always assume 

x n > and x n = 0. (30) 

Let B and B' be the Dn^ crystals. 

Definition 17. For each pair of variables x — (x\, . . . ,x\) G B,y = (j/i, . . . , y x ) G B', define the set 
of variables z^,z^ (0 < i < 2n - 2), yP,y?(l <i<n-l), x^x^y^y'^l < i < n) as follows. 



1. Set 



z 



(°)=y„, z^=y n . (31) 
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2. Define ^"-^""^t/P^P^ (1 < i < n - 1) as 



(32) 



7( ft, ^ n - 1 -*', z^ 1 -), y<) ari ) 

= ( z^\ yp, yp, *(»"*), y[ ), 

by descending recursion on z. Here the function 7 is given by Definition 16. 

3. Define z ^- l+i \ z^- 1+l \ x\,x\,y' r (1 < i < n - 1) as 

7 ( z<- 2 +*\ yp, yp, z("- 2+J ), ) 

= ( x[, ^"-i+O, ^, % ), 

by recursion on z. 

4. Set 

n-l 

< = z( 2 - 2 ), 34 = z( 2 - 2 \ y' n =£(x)-^(y'i+y'd, Vn = 0- (34) 

i=i 



(33) 



Remark 18. Let us consider the case when x G -Bi[n] and y G B\. Then the numbers represented 
by the variables z^ etc. in Definition 17 are equal to the numbers of items in the particle antiparticle 
description in Sec. IV B. More precisely these items appear within the time evolution process by T 
in Theorem 13 for the case N — 1. See Table III. In the table U and tj arc defined as follows: For 
1 < i < n— 1 we let U (respectively U) denote the time just after the z-carrier (respectively z-carrier) has 
passed, where the automaton state is given by /Q • • • K n -\Vy (respectively !€-■■■ /CyQ/Ci • • • K n -\Vy). 



Table III: The correspondence between the variables in Definition 17 and the items in the particle 
antiparticle description. 



Variables 


Items at time t% 


Variables 


Items at time U 


z (n-i) 


empty boxes 


~(n— 


empty boxes 


jin-i) 


boxes with bound states 


j(n-l+i) 


boxes with bound states 


yp 


boxes with z-balls 


A 


boxes with i- balls 


VP 


boxes with z-balls 




boxes with z-balls 


y'i 


balls in the z-carrier 




balls in the z-carrier 



C Proof 

We now give the proof of Theorem 13 that we have promised at the end of Sec. IV. It is obtained from 
the following. 

Theorem 19. Let x = (x\, . . . ,x{) G B,y = (j/i, . . . ,y x ) G B 1 be a pair of variables, and suppose the 
condition (30) on x. Let u,, Wi, v* , v^ 1 , v° n be the functions defined in Sec. V A, and let x\,x' i ,y' i ,y' i (\ < 
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i < n) be the variables given by Definition 17. Then the following relations hold 



x[ 




yi 


+ 


v a 


- v u 








x 'i 




Vi 


+ 


Vi-l 


- Vi+Wi-Wi-i 


(2<i< 


n — 


1), 


X n 




y n 


+ 


Vn-l 


- v° n 

u n > 








f'j 




y~i 


+ 


Vi-l 


- Vi 


(1 < i < 


n), 




y'i 




Xi 


+ 




- "i 


(1 < i < 


n), 




y'i 




Xi 


+ 




- «*, 








y'i 




Xi 


+ 


V*-l 


- v* +Wi- Wi-i 


(2<i< 


n — 


1), 


V'n 






+ 


<-l 


- v a " 









(35) 



A proof of Theorem 19 will be given after the following two lemmas. 

Lemma 20. Let A = y i7 B = min(x i+ i, y i+1 ), C = v* + mm(x i+ i,y i+1 ), D = yi and E = Xi in 
-y:{A,B,C,D,E)^{F,G,H,I,J). Then 

F = min(a; i ,y l ), 

G = Vi- m\n{x l ,y, i ) + min(a; i+ i, y i+1 ), 
H = -mm(x l ,y i ) + min(a; i+ i , y l+1 ) + y, + v* - v*_ x , 
I = v*_ 1 + min(a; i ,f l ), 
J = x i + v*_ 1 -v*, 

for 1 < i < n — 1 . 

In what follows we set wq = 2vq. Note that we have v* } = vq from Table II and w\ = vq + v^ 1 from 
Lemma 22 in the Appendix. 

Lemma 21. Let A = Vi-i + min(a; l , yA, B = y~i- min^,^) + min(x i+ i, y i+1 ), C = - min^fj + 
mm(x l+ i,y i+1 )+y i +v*-v*_ 1 ,D ^ min(x l ,y i )+v i -i+v*_ 1 -w l - 1 , and E = x~i inj: (A,B,C,D,E) ^ 
(F,G,H,I,J). Then 



F = 


Vi+Vi-i - Vi, 




G = 


v % + mm(x i+1 ,yi 




H = 


mm(x i+1 ,y i+1 ) -+ 


- Vi + v* - un, 


I = 


y { + vn - Wi-i - 


Vi + Vi-l, 


J = 


x l +w i - Wi-i - 


v* +v*_ 1 , 



for 1 < i < n — 1 . 

We shall give proofs of these lemmas in the Appendix. 

Proof of Theorem 19. Suppose i = n — 1 in Lemma 20. Then we have B — y n — z^ and C = 
y n = because of (25), (30), and (31). Then by comparing (32) with Lemma 20 we see that 
2^,1/^1,2/^-1,2^, and y' n _ 1 in (32) should be equal to F,G,H,L, and J. Thus the expression for 
y' n _i was obtained. The expressions for y[ (1 < i < n — 2), as well as those for z^ n ~ l \yP , yp , z^ n ~^ 
will be obtained by descending recursion on i, where one uses F and L as B and C in the next step. 
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Since B and C in Lemma 21 are equal to G and H in Lemma 20, we have B = yY and C = yf in 
Lemma 21. We also see that when i = 1 we have A = z'" -1 ) and I? = min(xi,y 1 ) = z'™ -1 ) in Lemma 
21, from the result obtained in the preceding paragraph. Then by comparing (33) with Lemma 21 we 
see that x[, z^ n \z^ n \ x[, and y\ in (33) are equal to F, G, H, I, and J in Lemma 21 if i = 1. Thus the 
expressions for x' l7 x' 1 ,y' 1 were obtained. The expressions for x' i ,x' i ,y' i (2 < i < n — 2) will be obtained 
by recursion on i, where one uses G and H as A and D in the next step. 

Then from (34) we can obtain the expressions for x' n and x' n in (35) since we have (30) and (25). 
It is clear that the relation y' n = = x n + v* l _ 1 — holds. Then the expression for y' n in (35) is 
obtained from the condition £(y') — £(x). The proof is completed. □ 

Finally we give the proof of Theorem 13. 

Proof of Theorem 13. If we impose the condition (30) on the defining relations (10) in Definition 5 then 
their right hand sides become those of (35) because of the existence of the limiting functions defined 
in Sec. V A. Then Theorem 19 tells that the numbers represented by x! i ,x' i ,y' i ,y' i in Definition 5 are 
equal to those by the same symbols in Definition 17 under this condition. Then according to Remark 
18 we see that the time evolution T by Definition 11 is identical to the T in Theorem 13 for the case 
N = 1. The assertion for the case N > 1 follows immediately by repeated use of this case, where we 
adopt the final states of the carriers for Bi i in (19) as their initial states for B li+1 . □ 

ACKNOWLEDGMENTS The author thanks Atsuo Kuniba and Akira Takcnouchi for a collabo- 
ration in the previous work. 

APPENDIX : PROOFS OF THE LEMMAS 

Before giving the proofs we present some relations between piecewise linear functions. They are used 
in the main text and in Appendix. These relations have been essentially obtained in Rcf. [13]. 

Ultradiscretization [18, 20, 21] is a procedure to derive an equation of piecewise linear functions 
from an equation of totally positive (i.e. having no minus sign) rational functions. It is realized as a 
transformation that replaces +, x and / by max(min), + and — , respectively. As the ultradiscretization 
of Lemma 4.12 of Ref. [13] we have the following. 

Lemma 22. The following relation holds 

W^Vo + Vf 1 . (Al) 

This lemma is used just above Lemma 21. The next lemma is obtained from the formulas (4. 23)*, (4. 23), 
and (4.24) of Rcf. [13] by the ultradiscretization. 

Lemma 23. The following relations hold 

max{V*J(x),e{x)+x i -y i }=max{xi-y i + V*_ 1 ,e(y)J(y) + Xi-y i }, (A 2) 

max{Vi,£(y),£(y) + y l - x,} = max{f 4 -x, + V l - 1 ,£(x),t(x) xj, (A 3) 

max{V n ,£(y) + X} = max{y„_ 1 +y n - z n _i -x n + V n - 2 ,£{x) + X}, (A 4) 

where 1 < i < n - 2 and X = y n -y n + (y n _ 1 + y n - x n -\ - x n ) + - 

Here we write (x)+ for max(x, 0). This lemma will be used in the proof of Lemma 14. 
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Proof of Lemma 10 

Proof. We derive the expression for Vj. First we suppose i ^ n — l,n. Then under the specialization 
(17) we have V t = maxi<_ 7 -< n _i{a i) j,a^} with 



£(x) for j < i, 
£(y) for j > i, 



a 'i,j = e ( x ) + ^2(Vk -x k )+Vj- 
fc=l 

Since £*»,.;(<») < o4i and a ij ( >i ) < aj ; „_i, we can drop off aj ; j(^„_i) in the max. Thus we obtain the 
desired result from (2) with i = l. Now we suppose i = n — 1 or n. Then we have 



auj = £{x) 



for i = n — 1, 

max(-y„, -x n -i) for i = n, 



fc=l 

Since ajj < a- j, we can drop off ctjj in the max and obtain the desired result. 

We derive the expression for V*. Note that if i = 0, n it has already been proved since V^* = for 
i = 0,n. Suppose i 7^ 0, n. Then under the specialization (17) we have V* — maxi<j< n _i{a| ., (c^ •)*} 
with 



+ i(#fc-J/k) fori < «, 
^(z) + Efe=i+i (j/fc - arfe) for i > 



( a i,j)* = S j,n-ix n + £(y) + ^2(x k - y k ). 



fe=i 



Since (c^ j(^ n _i))* < we ca n drop off (ofj^^^n-i))* in the max. The remaining candidates are 

a ti+i =t(x) +Vi+i, 

a i,i+2 = a i,i+l ~ x i+l + Vi+2-, 

1 

a i,n-l = a i,n-2 ~ x n-2 + Vn-1, 
( a i,n-lT = a i,n-l - x n-l + Vn, 
= «n-l)* ~Xn +V1, 
a i,2 = a i,l -X1+V2, 

Thus wc obtain the desired result. 
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We derive the expression for Vq 1 . Under the specialization (17) wc have V^ 1 = maxi<j< n _i{aQj, (ct' j) ai } 
with 

a 0j = %) +xi-{l- 5j.i)yi + 5j tn -ix n , 
j 

k=2 

Since £*oj(^n-i) — a o\ anc ^ ( a 'o.i) ai — ( a o.2)' Tl > we can drop off a^ 1 .^ and (ag.i)' 71 in the max 
and obtain the desired result. 

We derive the expression for V£ n . Under the specialization (17) we have V° n = maxi<j<„_i {a£ n . , {a' n j ) an } 
with 

a nJ = t( x ) + Vn ~ (1 - Sj, n -l)x n , 

= + y« - x„ + - x k ) + yj. 
k=i 

Since &n"j(^n-i) — ( a n.i) a ™ > we can drop off c^™-^^ in the max and obtain the desired result. 

The expression for W n _i is derived from (13), and that for W i ^ n _i~ ) is from (12) and the following 
lemma. □ 

Lemma 24. Let Pi be the function that was defined by (2). Then 

Pi+i > Pi — Hi- 
Proof. It is easy to see that 

P i -y i = max {A,}, P i+ i = max {xi + AA, 

l<j<n 2<j<n+l 

where 



Aj = ^(yt+k - x i+k -i). 



k=i 

The claim of the lemma holds since we have Xi > and Xi + A 2 = j/i+i > = A\. □ 
Proof of Lemma 14 

Proof. By definition we have V° n = maxi<j<„_i{a^™-, (a' n j) CT "}, where 
(anjY" = max(5 j>n _i(a;„ -y n ),y } - xj) 

n-l 

+ £(x)+ ]T {Vk-xk) + ]Jn X n , 

k=j+l 

( a 'n,jV n = max((5 J>l _i(f„ - x n ),Xj - yj) 

n-l j 



+ + ^2(Vk- x k) + ^2(vk - x k ) +y n - x n . 



fe=l fe=l 
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In the limit Km* the only element that survives in the max is (a nj „-i) CT ", which yields v^ n — y n — V n - 
In the same way the relation 

v n-2 = Vn-1 +Vn- »n-l + (x n -l ~Vn-\ + Vn) + , (A 5) 

can be obtained by a direct calculation. Then from (A 2) we see that the i>*'s exist and the relation (26) 
holds for 1 < i < n — 2 by descending induction on i. Since Vq = Vq (Table II) we have v = Vq. Then 
from (A 3) and (A 4) we see that v^s exist for i ^ n — 1, and that the relations (27) for 1 < i < n — 2 
and (28) hold, by induction on i. Since v n and v° n exist, we see by (16) that the function also 
exists and equals to max{v„ — y nl —y n }. Substituting (28) into v„-i = max{n„ — y n , —y n } we obtain 
(27) for i = n— 1. From * of (16) we obtain v^_ x = y n —y n - Then from (A 5) and max{— y n , —y n } = 
we obtain (26) for i = n — 1. From <j\ of (A 2) the existence of Vq 1 can be verified. The proof is 
completed. □ 

Proof of Lemma 15 

Proof. For 1 < i < n — 2 the relations follow immediately from (12). We consider the case i = n — 1. 
We obtain u>„_! = v n + from (13). First we suppose x„_i < y n _i- Then from Lemma 14 we 
obtain = y n _i + max{— x n -i + w„-2, —%n-i}- It yields the relations v n + v° n = v n -i + y n and 

RHS of (29) = w„_i + y n - y n + x n -i - min^-! - y n ,x n -i - v n - 2 - 

= v n -i +y n - min{0,f„ + (y„_ x - x n -i) + (z„_i - x n -i - v n - 2 )+} 
= v n -i + y n . 

Thus the assertion of the lemma was proved in this case. Now we suppose x„_i > y n _i- Then 

RHS of (29) = «„_! + y n - y n + y n _ x 

- min{x„_i - (x n _i - y n _ x - y n )+ -y n ,x n -i - v n - 2 - v n -i} 
= ma,x{(y n + y n _ 1 - z n _i)+ + v n -i,v n - 2 +V n -i -%n-i + Vn -y n }- 

From Lemma 14 we have v n -\ = (y n _i — x n -\ + v n - 2 ) + . Therefore 

RHS of (29) = max{(y„ + y n _ 1 - x n -i) + , 

y n _ l - x n -\ +y n + v n - 2 + max{-y„, -y n , y n _ 1 - x n _i}} . 

The last expression gives v n + v° n since the inner max vanishes. The proof is completed. □ 
Proof of Lemma 20 

Proof. The expressions for F and G are given by definition. The expression for / is given by I = 
C + D - H, and that for J is by J = I + B - C + E - F. Thus it suffices to prove H. We have 

H = min(C*, B + E-F) 

= min(a;j + i,37 i+1 ) +mm{v*,-min(xi,y i ) + Xi} 

= mm(x i+ i,y l+1 ) + v* - min(a; l ,f i ) + Xi - max {v* , - min(a; i , + x t } . 
Then by (26) we obtain the desired result. □ 
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Proof of Lemma 21 

Proof. The expression for G is given by G = A + B — F, that for I is by / = C + D — H, and that for 
J is by J = I + B — C + E — F. Thus it suffices to prove F and H. For F we have 

F = mm(A, E) 

= Vi +Vi-i + min{min(x J ,f i ) ~y i ,x l ~y i - Vi-i) 

= Vi+ Vi-i - max{-min(.T J ,yJ + y { , -x { + y i + Vi-i} . 

Then by (27) we obtain the desired result. For H we have 

H = min(C, B + E-F) 

= min(a; 4+ i,y i+1 ) - min(a; J ,y i ) + mm{v* - v*_ x +y l ,v l - +x,}. 

Then by (29) we obtain the desired result. □ 
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Figure 1: The loading and unloading process by the carrier of balls with index a. A gray ball represents 
the bound state of a particle and an antiparticle. It is assumed that b ^ a, a. 
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Figure 2: The meaning of the map 7 in Definition 16 by the particle antiparticle description of the 
automaton. 
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